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We present a magneto-infrared spectroscopy study on a newly identified three-dimensional (3D) Dirac 
semimetal ZrTej. We observe clear transitions between Landau levels and their further splitting under magnetic 
field. Both the sequence of transitions and their field dependence follow quantitatively the relation expected for 
3D massless Dirac fermions. The measurement also reveals an exceptionally low magnetic field needed to drive 
the compound into its quantum limit, demonstrating that ZrTe 5 is an extremely clean system and ideal platform 
for studying 3D Dirac fermions. The splitting of the Landau levels provides a direct and bulk spectroscopic 
evidence that a relatively weak magnetic field can produce a sizeable Zeeman effect on the 3D Dirac fermions, 
which lifts the spin degeneracy of Landau levels. Our analysis indicates that the compound evolves from a Dirac 
semimetal into a topological line-node semimetal under current magnetic field configuration. 
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3D topological Dirac/Weyl semimetals are new kinds of 
topological materials that possess linear band dispersion in 
the bulk along all three momentum directions[l-7]. Their 
low-energy quasiparticles are the condensed matter realiza¬ 
tion of Dirac and Weyl fermions in relativistic high energy 
physics[8, 9]. These materials are expected to host many 
unusual phenomena[10-12], in particular the chiral and ax¬ 
ial anomaly associated with Weyl fermions[3, 13-15]. It is 
well known that the Dirac nodes are protected by both time- 
reversal and space inversion symmetry. Since magnetic field 
breaks the time-reversal symmetry, a Dirac node may be split 
into a pair of Weyl nodes along the magnetic field direction 
in the momentum space [16-18] or transformed into line- 
nodes[17, 19]. Therefore, a Dirac semimetal can be consid¬ 
ered as a parent compound to realize other topological variant 
quantum states. However, past 3D Dirac semimetal materi¬ 
als (e.g. Cd 3 As 2 ) suffer from the problem of large residual 
carrier density which requires very high magnetic field (e.g. 
above 60 Tesla) to drive them to their quantum limit[20, 21]. 
This makes it extremely difficult to explore the transformation 
from Dirac to Weyl or line-node semimetals. Up to now, there 
are no direct evidences of such transformations. 

ZrTes appears to be a new topological 3D Dirac material 
that exhibits novel and interesting properties. The compound 
crystallizes in the layered orthorhombic crystal structure, with 
prismatic ZrTe f , chains running along the crystallographic a- 
axis and linked along the c-axis via zigzag chains of Te atoms 
to form two-dimensional (2D) layers. Those layers stack 
along the /?-axis. A recent ab initio calculation suggests that 
bulk ZrTes locates close to the phase boundary between weak 
and strong topological insulators [22]. However, more recent 
transport and ARPES experiments identify it to be a 3D Dirac 
semimetal with only one Dirac node at the T point [23]. In¬ 
terestingly, a chiral magnetic effect associated with the trans¬ 
formation from a Dirac to Weyl semimetal was observed on 


ZrTes through a magneto-transport measurement [23]. Our re¬ 
cent optical spectroscopy measurement at zero field revealed 
clearly a linear energy dependence of optical conductivity, be¬ 
ing another hallmark of 3D massless Dirac fermions [24]. 

In this letter, we present magneto-infrared spectroscopy 
study on ZrTes single crystals. We observe clear transitions 
between Landau levels and their further splitting under mag¬ 
netic field. Both the sequence of transitions and their field 
dependence follow quantitatively the relation expected for 
3D massless Dirac fermions. Furthermore, the measurement 
reveals an exceptionally low magnetic field (about 1 Tesla) 
needed to drive the compound into its quantum limit. Both 
facts demonstrate that ZrTes is an extremely clean system and 
ideal platform for studying 3D Dirac fermions. The presence 
of further splitting of Landau levels, which has never been 
observed in 2D massless Dirac fermions, e.g. graphene, pro¬ 
vides direct evidence for the lifting of spin degeneracy of Lan¬ 
dau levels, an effect being linked to the transformation from a 
Dirac semimetal to a line-node or Weyl semimetal. Our the¬ 
oretical analysis indicates that the former one is more likely 
realized in the present magnetic field configuration. 

Ligure 1 shows the reflectance spectra under different mag¬ 
netic field R(B) renormalized by the zero field reflectance /((()) 
in the far- and mid-infrared region. Lor the lowest magnetic 
field (1 T), a series of peaks could be clearly resolved, which 
keep growing more pronounced and shift to higher energies 
when the field strength B increases. In optical reflectance 
measurement, such peak features usually come from the in¬ 
terband transitions. Since these sharp peaks emerge in the re¬ 
flectivity only by applying magnetic field, it is natural to con¬ 
nect them to the Landau quantization of 3D Dirac electrons. 
Thus the peaks should stem from electronic transitions con¬ 
necting different Landau levels. Significantly, the first broad 
peak, which appears at the lowest energy, gradually split into 
four narrow peaks as B increases. This character is quite in- 
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FIG. 1: The relative reflectivity of ZrTes under magnetic field, as a 
function of energy. The spectrum are shifted upward by equal inter¬ 
val corresponding to different B. 


FIG. 2: The wave number dependent relative reflectivity R(B)/R( 0) 
under magnetic field of 2 T. n = 0... 5 point to six different peaks 
in sequence. The inset shows the linear dependence of the transition 
energies between Landau levels on ^fn + yn+ 1. The dash line is a 
guide to the eyes. 


triguing and has never been observed ever before, which will 
be explained in detail later. For the sake of convenience, they 
are marked by the numbers “1, 2, 3, 4” respectively at the top 
of Fig. 1. However, the splittings of other peaks located at 
higher energies are rather vague in this plot. 

In order to verify our speculation on the origin of the emerg¬ 
ing peaks and capture the underlying physics of 3D Dirac 
semimetal, we examine the sequence of peaks observed at low 
field. For a 3D system, the band structure would transform 
into a set of ID Landau levels by applying strong enough 
magnetic field, which are only dispersive along the field di¬ 
rection. Theoretical calculation on an isolated Weyl point 
has suggested that the magneto-optical conductivity is consti¬ 
tuted of a series of asymmetric peaks lying on top of a linear 
background [25, 26]. Especially, the peaks associated with 
allowed interband transitions in the Landau level structure oc¬ 
cur at o) oc \fn + V« + 1, corresponding to transition from 
L_„ to L n+ 1 or from L_(„+p to L„, where L n represents for 
the nth Landau level. This conclusion applies to 3D massless 
Dirac fermion as well, because massless Dirac fermion can be 
thought as two sets of Weyl fermions with opposite chirality. 

A linear rising optical conductivity has been revealed in our 
previous zero field spectroscopic experiment on ZrTes single 
crystal [24] which already provides a strong evidence for 3D 
massless Dirac or Weyl femions. We blew up the results of 
B — 2 T as displayed in Fig. 2, in which six peaks could be 
clearly resolved. The positions of these peaks are identified 
to be about 202, 480, 628, 748, 856, 937 cnG 1 in sequence. 
The energy ratios of the peaks observed here can be approx¬ 
imately reduce to 1 : 1 + V2 : V2+ V3 : V3+ V4 : V4+ V5 
: V5+ V6, in nearly perfect accordance with the predicted 
massless Dirac semimetal behavior. This results in a linear 
dependence of the transition energy between Landau levels 
on yfn + V« + 1, as shown in the inset of the figure. From this 
observation, the first peak can be unambiguously determined 


to correspond to n— 0, ascribed to transitions from Lq to L\ and 
L_i to L{). Only when the chemical potential lies in between 
L\ and L \, can this transition be clearly visible. [26] This is 
quite exciting because it demonstrates that the quantum limit 
could be easily approached by magnetic field as low as 1 T, 
where the n - 0 peak is distinctively observed. As a con¬ 
trast, the quantum limit of the well-known Dirac semimetal 
Cd 3 As 2 can not be reached with magnetic field lower than 65 
T [20, 21]. It indicates that ZrTes compound is extremely 
close to an ideal Dirac semimetal, with the chemical poten¬ 
tial lying in the vicinity of Dirac point, and meanwhile it is 
an extraordinarily clean system. We anticipate that our find¬ 
ing of easy access of quantum limit will motivate many other 
experimental probes on this compound. In Supplemental In¬ 
formation [27] we perform more detailed analysis of this peak 
sequence and obtain estimates of the average a (--plane Fermi 
velocity v_i_ = y/v a v c ~ 4.84 x 10 5 (m/s) and a vanishingly 
small Dirac mass \m\ ~ 2 (cm -1 ). 

In Fig. 1, where R(B)/R( 0) was shifted by equal interval 
with regard to increasing magnetic field, it is noted that the 
peak positions evolve in a way much alike the parabolic fash¬ 
ion as B. To further illustrate the characteristic features of the 
Landau levels, we plot R(B)/R( 0) in a pseudo-color photo¬ 
graph as a function of Vb. It is clearly seen in Fig. 3 that the 
wave numbers of the peaks are basically linear proportional 
to Vb. The dashed lines are instructions for eyes, whose in¬ 
tercepts at 0 T are all absolute zero. For a single massless 3D 
Dirac node, the nth Landau levels caused by magnetic field 
are dispersive only along the field direction (see Supplemen¬ 
tary Information[27] for more details), with doubly degener¬ 
ate n + 0 levels E n (k\\) = sgn(n) ■ ^2 v\eBti ■ \n\ + ft 2 vjk^, 
and Eo(k\\) = ±hv\\k\\ for n — 0. If neglecting the disper¬ 
sion along the magnetic field direction, then E„ oc Vb. As 
a comparison, for free electron systems, the magnetic induced 
Landau level obeys E n — (n + j)ha> c , where a> c is the the 




































FIG. 3: The pseudo-colors photograph of the relative reflectivity 
R(B)/R( 0) as functions of wave number and ^^B. The dashed lines 
are linear fittings of the peak energies dependent on V5. 


cyclotron angular velocity and proportional to B instead of 
Vs. Additionally, the Landau level energy of massive 3D 
Dirac fermions in topological insulator Bi 2 Se 3 is reported to 
be in linear scale with B as well[28]. Therefore, the \fH de¬ 
pendence in Fig. 3 intensively implies again the characteristic 
property of 3D massless Dirac fermions in ZrTes under mag¬ 
netic field. 

In addition to the sequence of peaks at low field, the peak 
splitting shown in Fig. 1 could also be well resolved in Fig. 3. 
The n — 0 peak evolves into four peaks at high magnetic field, 
with two very pronounced ones at lower energies and two rel¬ 
atively weak ones at higher energies. They are labeled as ”1, 
2, 3, 4” respectively in accordance with Fig. 1 The splittings 
of the rest peaks are, although too vague to be precisely iden¬ 
tified, but for certain to exist. The n = 1 peak, arising from 
transitions from L 2 to L\ and L_i to Li, seems to split into 3 
peaks. Such splitting has never been observed in 2D massless 
Dirac fermion system, for example graphene[29]. 

We now explore the underlying mechanism for the splitting. 
We will show that this splitting can be naturally explained by 
the Zeeman effects of magnetic field on 3D Dirac fermions. 
The effect of Zeeman field on Dirac semimetals has been thor¬ 
oughly studied by Burkov el al. [17], It was pointed out that 
Zeeman field may split the Dirac node into two Weyl nodes, 
or transform the Dirac node into “line-nodes”[17]. Further in¬ 
cluding orbital effects of magnetic field will generate Zeeman- 
split Landau levels, schematically shown in Fig. 4. Here we 
will briefly explain the two possible scenarios and leave the 
details in the Supplementary Information [27]. 

The first scenario we consider is the “line-nodes” picture, 
with the resulting Landau level structure depicted in Fig. 4(a). 
The doubly degenerate Landau levels E n (k\\) for n 4- 0 from 
3D Dirac fermions will be split into E n ± ~ E„(k\\) ± g^BB/ 2, 
where g is the average Lande g-factor for the conduction and 



FIG. 4: Two proposed scenarios. The left panel (a) shows the split¬ 
ting of Landau levels simply caused by Zeeman effect. The strong 
spin-orbit coupling could lead to a mixing of spin up and down com¬ 
ponents. So the split ”+” (solid lines) and (dashed lines) levels do 
not have pure spin up and spin down components. This would allow 
inter Landau level transitions to occur with opposite signs, however, 
weak in intensity. The purple dashed line is the Fermi energy E F . 
Zeeman field without orbital effect will transform 3D Dirac node into 
“line-nodes” in this case[17]. The right panel (b) represents for the 
Landau levels of a Weyl semimetal induced by the magnetic field. 
The crossing points of n = 0 Landau levels with the horizontal axis 
are the would-be Weyl points, if the orbital effect of magnetic field 
is ignored. The solid lines and dashed lines represent the two sets of 
Landau levels from the two Weyl nodes with opposite chirality. Tran¬ 
sitions between Landau levels of opposite chirality will have weaker 
intensity. In both panels, the red solid arrows represent for transitions 
between Landau levels of the same spin/chirality, whereas the green 
dashed ones indicate that of different spin/chirality. 


valence bands of 3D Dirac fermion. The n — 0 Landau levels 
£o(^||) = ±7iV||fc|| will mix around k\\ = 0 and open a gap of 

size gj-igB there, and become Eq ± ~ ± ^h 2 v 2 k 2 + ( gfisB/2 ) 2 . 
The split Landau levels are labeled by “spin” indices '+’ and 

’ in Fig. 4(a), which indicate that the states are of spin up 
or spin down at k\\ = 0. However with k\\ + 0 the split Landau 
levels are not purely spin up or spin down due to strong spin- 
orbit coupling. Optical transitions between levels of the same 
“spin” indices can happen at k\\ = 0, and produce strong peaks 
in optical conductivity (thus reflectivity) [26], but transitions 
between levels of opposite spin indices will be suppressed at 
k\\ = 0, leading to weak and broad peaks. Therefore the origi¬ 
nal n — 0 peak will split into two strong peaks “1” from Lo- 
to L| _, and “2” from L_i,+ to L (l f , and two weak peaks “3” 
from Lo _ to L\ + , and “4” from L_j _ to Lo,+ - The peak “1” and 
peak “2” can have different energy if the chemical potential is 
not in the gap between Lq _ and Lo,+ as depicted in Fig. 4(a), 
or if the conduction and valence bands of 3D Dirac fermion 
have different ^-factor (see Supplementary Information [27]). 
In any case the splitting between peak “3” and peak “1”, and 
between peak “4” and peak “2” will be about giisB. Based on 
previous ab initio results[22] on ZrTes and the experimental 
setup we conclude that this scenario is the most likely expla¬ 
nation of our observation (see Supplementary Information[27] 
for details). 

The second and more interesting scenario is the “Weyl 
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FIG. 5: The magnetic field dependent energy difference of EyE\ 
and £ 4 -£ 2 - They represent the Landau level splittings for the con¬ 
duction and valence bands, respectively. 

nodes” picture illustrated in Fig. 4(b). In this case the Zee- 
man field effectively shifts the wave vector parallel to field 
by ±g/jBB/2hv\\, where the ± sign depends on the chirality 
of Weyl fermions. The degenerate n f 0 Landau levels be¬ 
come E n ± (k\\) ~ E„(k\\ + g^igB and the n = 0 Landau 
levels become £o,±(k||) ~ - g^sB/ 2. In this scenario 

transitions between Landau levels of the same (different) chi¬ 
rality will have strong (weak) intensity. The original n - 0 
peak will also split into two strong peaks “ 1 ” from Lq.-s to 
L\ - s , and “2” from £_i _,j to Lo.-r, and two weak peaks “3” 
from Lo.-s to L\ s , and “4” from £ i jS to Lq.-j, where s - +. 
The peak “1” and peak “2” can have different energy if the 
chemical potential is not at the charge neutrality as depicted 
in Fig. 4(b). However the splitting between peak “3” and “1” 
will in general not be linear in B, unless the conduction and 
valence bands have very different g-factors. According to our 
analysis[27] this scenario is more likely to happen when the 
field is applied along crystal c-direction. 

From the peak positions of “1”, “2”, “3” and “4”, we can 
immediately obtain the dependence of the split energy as a 
function of magnetic field. Figure 5 displays the magnetic 
field dependence of £4 - £3 above 6 T and EyE \ above 13 T, 
respectively. The energy positions for peak ”3” could not be 
well resolved below 12 T, so the energy difference of £3 - E\ 
is plotted only at high magnetic field. Obviously, both £4 - £2 
and £ 3-/41 exhibit good linear dependence, in better agree¬ 
ment with the first “line-nodes” scenario. The energy split¬ 
ting is roughly 10.5 cm '/T (~1.3 meV/T) for £4 - £2 and 
7.4 cm 1 /T (~0.92 meV/T) for £3 - E\. According to previ¬ 
ous discussions, this leads to estimates of the average factor 
being 22.5 or 15.8, from the £4 - £2 or £3 - E\ splittings 
respectively. 

It was known that the g-factor reaches about 37 for 
Cd3As2[30], which is even bigger than our estimated ^-factor 
for the ZrTes compound. We do not yet have a good explana¬ 
tion for the discrepancy between £ 4-£2 and £3 - £ 1 splittings. 
Based on our analysis in Supplementary Information[27], this 
could be due to the broad nature of the weaker peaks, and the 


current analysis may overestimate the ^-factors for this rea¬ 
son. 

The transitions associated with £_2 —> L l and £ 1 —> £ 2 are 
much more complex. Analogy to the transitions between £ () 
and £ ± i, the selection rule permissive ones (indicated by red 
solid arrows) should be more pronounced than those between 
opposite spin orientations (the dashed green ones). Consider¬ 
ing possible different splittings of the valence and conduction 
bands, the n — 1 peak appeared in R(B)/R{ 0) is supposed to 
contains at least three components. 

In summary, by performing magneto-optical measurement 
on the single crystalline 3D massless Dirac semimetal ZrTes, 
we have clearly observed the magnetic field induced Lan¬ 
dau levels, evidenced by regular organized peaks shown in 
the renormalized reflectivity R(B)/R( 0). Particularly, the first 
peak is identified to be originated from the transitions between 
the zeroth and first Landau Levels, which reveals the Fermi 
energy lies very close to the Dirac point. The appearance 
of the first peak under magnetic field as low as 1 T demon¬ 
strates an exceptionally low quantum limit of ZrTes compared 
to other 3D Dirac semimetals, which provided an elegant plat¬ 
form to explore more intriguing non-trivial quantum phenom¬ 
ena. Of most importance, the fourfold splitting of the first 
peak yield direct and clear evidence for the release of spin 
degeneracy of Landau level, hence the transformation from a 
Dirac semimetal into line-node or Weyl semimetal. Our the¬ 
oretical modeling and analysis indicate that the former one is 
more likely realized in the present magnetic field configura¬ 
tion. 
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